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Abstract 

The objective of this research is to investigate the solder joint 
reliability of board-level 3D pyramid package-on-package 
(POP) stacked-die packages in transient drop analysis. Three 
kinds of lead-free solder alloy, Sn4Ag0.5Cu, Sn3.5AgO.5Cu 
and Sn3.5Ag are introduced to demonstrate the transient 
dynamic response subjected to JEDEC sine wave pulse- 
controlled drop test. In order to evaluate the structure of the 
interested area a strip model sliced from, the full test vehicle 
is introduced in this research. In addition, the submodel 
region is particularly chosen with the strip model by 
performing the cut boundary interpolation. The envelope of 
equivalent stress for the outermost solder joint off the end of 
the strip model is drawn to illustrate the potential failure 
sites and mechanism. The cut boundary of submodel is 
verified and the mesh density of submodel is examined. For 
a refinery mesh of submodel, parametric studies for 
structure and material are carried out to investigate the 
reliability of the outermost solder joint, and the results are 
summarized as design rules for the development of 3D IC 
packages. A series of comprehensive parametric studies are 
conducted in this research. 
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Introduction 

The tendency of advanced integrated circuit (IC) 
package techniques is characterized by lightness, 
thinness, small size and high performance. Hence, 
advance in the design and manufacture of the 
packaging technology need to keep pace with advance 
in the miniaturization of the electronic systems. The IC 
packaging technology plays an important role in 
determining the overall performance in the electronic 
system. The C4 (Controlled Collapse Chip Connection) 
packaging technique has been improved to ball grid 



array (BGA), flip chip (FC) and package-on-package 
(POP) technology in recent years [1-15]. In compliance 
with manufacturing process, the three dimensional 
(3D) packaging technique [1-2,4,7-9] presents the 
shortest leads with the minimized inductance as well 
as maximized frequency, and will become the leading 
packaging technologies in next generationO There 
have been many 3D stack-die POP packaging 
structures reported [2-3] and developed for their easily 
compliance with the current manufacturing process. 

The function of solder joints in the presented 3D POP 
structure are: electric interconnection, structural 
support and thermal conduction. In addition, the 
eutectic solder has been replaced by lead-free solder 
alloy in accordance with RoHS and WEEE 
environmental regulations. However, the primary 
difficulty encountered in these structures is the lead- 
free solder joints reliability (SJR) which is the 
probability of solder joints to perform their intended 
mechanical and electrical specifications over a given 
period of time and under a specified set of operating 
conditions. Four mechanisms of solder joint failure 
have been reported [13], such as tensile fracture due to 
mechanical overstress, creep failure caused by a 
permanent sustained load, fatigue resulted from a 
cyclic loading, and electronic migration owing to 
current crowding and working- temperature. Therefore, 
SJR studies, such as cycling temperature loading, 
electricity current crowding and transient drop effects, 
must take these failure mechanisms into consideration. 
SJR studies cover module-level and board-level solder 
joint reliability. Module-level SJR copes with the 
reliability of solder bumps within the package before 
board mounting. Board-level SJR deals with the 
reliability of solder joints on a package after board 
mounting. Board-level SJR for 3D POP package is 
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much more comprehensive because it comprises both 
the solder-to-package and solder-to-board interfaces. 

Previous works [2-3] demonstrated that two categories 
of board-level drop tests, namely free drop test and 
pulse-controlled drop test, were developed to 
investigate the transient dynamic structural behavior 
of solder joints. In order to evaluate the solder joint 
reliability of on-board various electronic packages, the 
Joint Electronic Devices Engineering Council (JEDEC) 
proposed the pulse-controlled drop test regulations, 
JESD22-B104, JESD22-B110 and JESD22-B111, 
respectively. For packaging body size less than 15mm 
x 15mm, JESD22-B111 regulation of the half sine pulse 
of 1500G acceleration within 0.5ms duration time is 
the commonly used test standard in the field. Since 
board-level reliability testing is more difficult to 
implement, an alternative scheme to analyse solder 
joint strengths and weaknesses is finite element 
analysis (FEA) using computer models. Few 
methodologies, such as the response spectra analysis 
[7], the input G method [1] and the support excitation 
scheme [4-5], were developed to investigate the pulse- 
controlled board-level drop test effects. Many recent 
works [2-9] have been reported that the support 
excitation scheme produces no rigid-body motion to 
the test vehicle and dramatically reduces the 
computational cost. 

The preliminary results show that the peak effective 
stress is located in the outermost solder joint off the 
end of strip model, which is the interesting region. 
Compared with the conventional numerical approach 
for drop tests, submodel technique provides both 
accurate and computationally economical solutions. 
Hence, the purpose of this research mainly 
concentrates on examining the transient dynamic 
response in the outermost solder ball for board-level 
drop tests based on the support excitation scheme. 

Constitutive Model 

Drop Test 

Previous researches revealed drop tests were 
conducted by using KD-DP1200 drop apparatus. FIG. 
1 demonstrates the drop test apparatus for the pulse- 
controlled board-level drop test based on JEDEC 
standard JESD22-B111. Detailed view of the drop 
system is schematic illustrated in FIG. 2, the drop table 
is freely released and along side guides from a certain 
altitude to hit the strike surface on the rigid base. The 
board-level test board is supported on the base plate 



with 4-corner fixed standoffs. Test board consists of 
fifteen representative samples bearing the solder joint 
of interesting POP stack-die packages. The base plate 
is welded on the drop table and the accelerometer is 
firmly attached to the base plate to monitor the peak 
acceleration and duration of half-sine impact pulse. 




t Strike Surface 
^ Rigid Base 



FIG. 1 SCHEMATICALLY ILLUSTRATION OF DROP TEST 
APPARATUS KD-DP1200 



Sample 
Standoffs 
Base Plate 



Drop TaMe 




FIG. 2 DETAILED VIEW OF THE DROP SYSTEM 

Support Excitation Scheme 

Based on the support excitation scheme[4-5, 7-9], the 
drop test system presented in FIG. 2 can then be 
transferred into a multiple degree-of-freedom (MDOF) 
of mass-damper-spring time dependent vibration 
system. It should be noted that the origin of the 
coordinate system is located on the top surface of base 
plate where the accelerometer is attached. FIG. 3 
depicts the corresponding MDOF time dependent 
vibration system. 
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FIG. 3 MDOF OF MASS-DAMPER-SPRING TIME DEPENDENT 
VIBRATION SYSTEM FOR DROP SYSTEM IN FIG.2 

In term of the force equilibrium along the X-direction 
on the MDOF system, equation of motion can be 
derived as follows 
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where m; is the mass, Ui is the displacement, dj is the 
damping coefficient, kj is the spring constant and F( t) 
is the impact force. The initial conditions are 

(2a) 



U X (0) = U 2 (0) = U 3 (0) = U 4 (0) = U 5 (0) = 
t/!(0) = C/ 2 (0) = C/ 3 (0) = t> 4 (0) = t/ 5 (0) = V 



(2b) 



where V is the initial impact velocity. After the drop 
table strikes the base table, the drop table, base plate 
and standoff are immovable. The aftershock 
displacements for the test board and sample packages 
can be observed from the base plate where the 
accelerometer is mounted. Another coordinate system, 
x is selected to observe the vibration of mi, mi and mi. 
The origin of x system is located at the observation 
surface and w; is the relative displacements. All the 
displacements can then be defined as 



Ui=U 3 + Ui 



(3) 



As mentioned earlier, Ui is the local coordinate system 
located on the top surface of base plate where the 
accelerometer is mounted. Notice that gravity is along 
the negative x-direction. Rewrite Equation (1) and the 
full equations of the system in x-direction can be 
expressed as 
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Rearrange Equation (4) to 
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However, the aftershock relative displacements 
remains only between mi and mi, and between mi and 
nt3, respectively. The board level drop test system can 
then be separated as two independent structure and 
Equation (5) can then be rewritten in two independent 
equations as 
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The initial conditions for drop test board with sample 
packages in Equation (6a) are 

M 1 (0) = M 2 (0) = M 1 (0) = li 2 (0) (7a) 
The initial conditions for drop table and base plate in 
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Equation (6b) are 

m 4 (0) = u 5 (0) = ii 4 (0) = u 5 (0) (7b) 
In FIG. 3, the acceleration is defined as 



G(t) = G m sin(— ) 



(8) 



The velocity and displacement are derived as follow 

r ct 7tT T TTt 

V(t) = )G(t)dt = LG m sm(—)dT = G m — (1- cos(— )) (9) 
"7 n 1 



D(t) = \V(t)dt = fG m — (l-cos(— ))dz 

U K T 



:G m— ( f Sln (— )) 



(10) 



In this research, Gm is 1500G and T is 0.0005s. 
Equation (6a) - (10) can be easily applied to the finite 
element model of board level drop impact test. 

Finite Element Model 

Finite element modeling based on ANSYS for a board- 
level test vehicle subjected to JEDEC drop test 
standard JESD22-B110 is revealed in this section. A 
strip finite element model sliced from the entire 
system is applied to satisfy the support excitation 
scheme. In order to save the computational time and 
remain the accuracy, submodelling technique around 
the interested solder joint is modeled. 

Strip Model 

Fifteen of 3D pyramid stack die POP package were 
mounted on the 132 mm x 77mm test board shown in 
FIG. 4 to study solder joint reliability on the dynamic 
drop effects in this study. Preliminary results show the 
maximum bending deflection which is found in the 
region near sample Us. This would imply the 
maximum amplitude occurred in this area and the 
peak effective stress on solder joint is reported. Thus, a 
sliced strip model is selected and applied to evaluate 
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the dynamic impact effects. FIG. 5 shows the finite 
element meshed strip model with boundary 
conditions for drop impact evaluation. FIG. 6 displays 
finite element 2D meshed strip model. Detailed 
dimensions of POP package used in this study are 
listed in TABLE 1. 



(Left) Fixed in x axe 



(Back) Fixed in z-axe 




CFront) Fixed in z-axe 



(Eight) Fixed in x, y, z axe ' 



FIG. 5 FINITE ELEMENT STRIP MODEL WITH BOUNDARY 
CONDITION FOR DROP TEST 



Epoxy Molding Compound 

BI Substrate 




FIG. 6 FINITE ELEMENT MESHED MODEL (2D FRONT VIEW) 

TABLE 1 DIMENSIONS OF POP PACKAGE USED IN THIS 
STUDY 



PoP Package 


Dimensions 


Die Size 


6 mm x 6mm 


Die Thickness 


4 mil 


Package Size 


15 mm x 15mm 


Die Attach Thickness 


1 mil 


Bottom Compound Thickness 


0.35 mm 


Upper Compound Thickness 


0.45 mm 


Bottom Substrate Thickness 


0.26 mm 


Upper Substrate Thickness 


0.21 mm 



FIG. 4 STRIP MODEL SLICED FROM THE ENTIRE TEST BOARD 



(note: 1 mil=0.0254 mm) 

Material Property 

TABLE 2 and TABLE 3 show the material properties 
for selected SnAgCu (SAC) alloys as well as polymeric 
materials used in this study. It should be noted that E 
is the Young's modulus, v is the Poisson's ratio, P is 
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the density and G is the shear modulus. 

TABLE 2 MATERIAL PROPERTIES USED IN THIS STUDY 



Material 


E(GPa) 


V 




Substrate/Test 
Board 


Transversely 
Isotropic(Refer to Table 3) 


1.91 


Sn4.0Ag0.5Cu 


43 


0.3 


8.42 


Sn3.5AgO.5Cu 


52 


0.36 


7.44 


Sn3.5Ag 


46 


0.3 


7.39 


Die 


131 


0.23 


2.33 


Die Attach 


2.7 


0.35 


3.15 


Compound 


28 


0.35 


1.89 



TABLE 3 TRANSVERSELY ISOTROPIC MATERIAL PROPERTY OF 
SUBSTRATE AND TEST BOARD 



Ex, Ey 

(GPa) 


Ez 
(GPa) 


GxZf Gyz 

(GPa) 


Gxy 

(GPa) 


Vxz / Vyz 


Vxy 


16.8 


7.4 


7.59 


3.31 


0.39 


0.11 



All the components are assumed linearly elastic except 
the solder joints. The material behavior of solder joints 
should reflect its large impact deformation at very 
high strain rate. Elastic-plastic tensile mechanical 
properties for selected SAC alloys were measured and 
the true stress-strain curves were illustrated in FIG. 7. 
Multilinear kinematic hardening (MKIN) material 
model is taken into account the effects of time- 
independent plasticity behavior when the solder ball 
suffers bending and dynamic phenomena. 
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FIG. 7 ELASTIC-PLASTIC TENSILE MECHANICAL BEHAVIOR 
FOR SAC ALLOYS 

Damping coefficients in Equation (6) can be formed by 
Rayleigh damping matrix 

[c]=a[m] +/?[&] (11) 

where [m] is the mass matrix and [k] is the stiffness 
matrix. The coefficients of a and are set to be 0.1 
and 0.1 in ANSYS. 



Submodel Technique 

Based on St. Venant's principle, the cut-boundary 
interpolation displacement method from coarse model 
is applied to obtain the results for the small interesting 
area. FIG. 8 presents the sliced strip coarse model and 
submodel of outermost solder ball off the end. FIG. 9 
illustrates the enlarged view of outermost solder ball 
in the submodel. In order to examine the potential 
fracturing area, the top surface of solder ball refers to 
the package side, while the bottom surface is referring 
to the test board side. 




FIG. 9 ENLARGED VIEW OF MAPPED MESH OF SOLDER BALL 



Results and Discussion 

Finite Element Analysis 

A general purpose software ANSYS based on finite 
element method is utilized to predict the dynamic 
drop impact reliability of SnAgCu alloy solder balls. 
Sn4Ag0.5Cu alloy was first applied to the material 
model to study mesh density and the stressed area on 
the solder ball. The predicted peak von Mise's stressed 
area was found in the lower level, the outermost 
solder ball where plastic failure was reported. FIG. 10 
demonstrates the enveloped von Mise's stress 
contours on the solder balls from center to the 
outermost along the strip model. 
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FIG. 10 ENVELOPED VON MISE'S STRESS CONTOURS ON THE 
SOLDER BALLS 

FIG. 11 illustrates the change of vibration amplitude at 
the maximum stressed observation point. Obviously 
the hilltop of the first amplitude is around 1.87 mm 
(downward) at 0.8 ms aftershock. In the second 
amplitude, the hilltop is around 1.23 mm (upward) at 
2.25 ms. However, the oscillation tends to return to the 
position of equilibrium, i.e. amplitude becomes zero 
due to the effect of Rayleigh damping. FIG. 12 presents 
the response von Mise's stress without time- 
independent plasticity behavior at the outermost 
solder ball. As it can be seen, von Mise's stress exceeds 
the elastic region in the first oscillation. Thus, the effect 
of time-independent plasticity behavior needs to be 
included in the material model. 




FIG. 11 PREDICTED AMPLITUDE VS. TIME AT THE 
OUTERMOST SOLDER BALL 




Submodel Technique 

Since the peak effective stress is located at the outmost 
solder ball in the lower level, it is of interest to 
evaluate the location of potential fracture in the 
outermost solder ball. The enveloped equivalent stress 
contours for coarse model (strip model) and submodel 
(outermost solder joint) were observed. The peak 
equivalent stress is found to be 79.4 Mpa in sliced 
coarse model, 74.22 MPa (688 elements) and 72.218 
MPa (17424 elements) in submodel. Comparison of 
coarse model and submodel is listed in TABLE 4. 

TABLE 4 COMPARISON OF COARSE MODEL AND SUBMODEL. 



Finite Element 
Model 


CPU Time 
(hr) 


Effective 
Stress (MPa) 


HD Storage 
(GB) 


Coarse model 
(162120 elements) 


44.3 


79.4 


45.8 


Submodel (688 
elements) 


1.4 


74.22 


9.8 


Submodel (17424 
elements) 


11.2 


72.218 


17.6 



FIG. 13 demonstrates the predicted von Mise's stress 
contours on the region near substrate side (top) and 
the peak von Mise's stress is 63.394 MPa. However, 
predicted peak von Mise's stress is 74.22 MPa in FIG. 
14, which implies that the potential plastic failure 
would occur near test board side (bottom). 
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FIG. 13 PREDICTED VON MISE'S STRESS CONTOURS ON THE 
REGION NEAR SUBSTRATE SIDE (TOP) 
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FIG. 12 PREDICTED VON MISE'S STRESS AT THE OUTERMOST 
SOLDER BALL 



FIG. 14 PREDICTED VON MISE'S STRESS CONTOURS ON THE 
REGION NEAR TEST BOARD SIDE (BOTTOM) 
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Lead-free Material 

Lead-free SnAgCu solder alloys have been widely 
used in the IC package industry in recent years. In this 
study, Sn4Ag0.5Cu, Sn3.5AgO.5Cu and Sn3.5Ag alloys 
are introduced to study the effect of lead-free solder 
on drop impact response. TABLE 5 lists the predicted 
peak amplitude and the predicted maximum effective 
stress during drop test. For Sn3.5AgO.5Cu and 
Sn3.5Ag, solder ball without copper is much more 
flexible. For Sn4Ag0.5Cu and Sn3.5AgO.5Cu, an 
increase in the weight percentage of silver (Ag) in the 
solder ball will result in a slightly increase in the 
amplitude. 



TABLE 5 PREDICTED AFTERSHOCK RESPONSE FOR DIFFERENT 
LEAD-FREE SOLDER ALLOYS 



Material 


Sn4Ag0.5Cu 


Sn3.5AgO.5Cu 


Sn3.5Ag 


Predicted 
Amplitude (mm) 


1.87 


1.84 


2.03 


Predicted Effective 
Strain (%) 


0.189 


0.185 


0.24 


Predicted Effective 
Stress (MPa) 


76.4 


65 


61.4 



Conclusions 

This paper demonstrates an insight into the plastic 
failure mechanism associated with dynamic drop 
impact analysis. Finite element predictions reveal the 
potential fracture region at the outermost solder joint 
in the lower level off the end of the strip model. Based 
on parametric studies in this study, some conclusions 
have been summarized as follows: 

1. The envelope equivalent stress profile has been 
reached maximum at the outermost solder joint in the 
lower level in the strip model. 

2. The peak equivalent stress, located at the surface 
near the test board side, implies that the solder 
fracturing initially occurs at this region. 

3. The reliability of Sn4.0Ag0.5Cu and Sn3.5AgO.5Cu 
alloys show the peak effective stress would be 
decreased as the ingredient of Ag decreases in the 
solder, which implies solder joint can absorb more 
heavy impact distortion energy as tin (Sn) is increased. 

4. For Sn3.5AgO.5Cu and Sn3.5Ag alloys, the peak 
effective stress increases as the ingredient of copper 
(Cu) increases in the solder, i.e. SnAgCu alloys has 
better reliability. 

5. For structure analysis, an increasing in the thickness 



of molding compound or a decreasing in the height of 
test board will produce a lower peak effective stress 
and improve the reliability of solder joints. 

6. An increase in the height of solder ball results in a 
larger deformation on the solder ball which is 
threatening the structure of solder ball. 
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